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Abstract 

The zero-point energy of a massless fermion field in the interior of two parallel plates in a D- 
dimensional space-time at zero temperature is calculated. In order to regularize the model, a 
mix between dimensional and zeta-function regularization procedure is used and it is found that 
the regularized zero-point energy density is finite for any number of space-time dimensions. We 
present a general expression for the Casimir energy for the fermionic field in such a situation. 
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Quantum Chromodynamics (QCD) is a non-abelian gauge field theory that describes the quark 
dynamics. The non-observation of isolated quarks, i.e., the fact that only colorless states are 
physically realized, introduced the concept of confinement. The mechanism by which this occurs 
is still unknown. The bag model is an attempt to conciliate the success of the quark model with 
confinement and its key point is the approximation that the bag is a classical spherical cavity 
with quarks and gluons confined. In a first approximation, the quarks and gluons move freely 
inside the bag, but are absolutely confined to its interior. Since any quantum field submitted to 
classical boundary conditions gives rise to the Casimir effect , we expect that confinement will 
give rise to the Casimir energy of the gluon and quark fields 0]. Over the past twenty five years 
many authors calculated the renormalized zero-point energy of different quantum fields in the bag 
model. The renormalized vacuum energy due to the scalar field in the bag model was presented by 
Bender and Hays |0 and more recently by Romeo and also Bordag et al. Q. The calculation 
for the case of the gluon field was presented by Milton [|| and also Romeo and for the fermion 
field Milton p| and also Bender and Hays obtained the Casimir energy in the bag model. 

In the renormalization problem, where we have to extract physically finite answers, two different 
physical situations must be studied. The first one is the case where the surface whereon the fields 
satisfy some boundary conditions divides the manifold into two parts, and the fields are present 
in both regions. The second one - the bag model case - is the case where the fields are confined 
to the interior of some region, i.e., there is no field outside 0. In this second case, since there are 
no exterior modes, if we use a cut-off or the Green's function method to regularize the zero-point 
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energy, divergencies appear and we will be able to absorb them by means of a renormalization 
procedure only if we introduce ad hoc contact terms. There are many attempts to solve this 
problem, but to our knowledge, this question is still open in the literature. Using the Weyl theorem 

it it not difficult to show that the asymptotic eigenfrequency distribution of any quantum field 
confined in a finite volume will present divergent terms proportional to geometric parameters of 
the region to which the field is confined. To keep in mind the problem, let us consider it for the 
case of the a scalar field. (Units in which h = c = 1 will be used throughout this letter. The 
metric tensor is taken to be r^^jy = diag{l, —1, —1).) 

Let us use the mode sum energy {E)^°^'^ as the total renormalized energy defined by: 

Jo 2 

where is the zero-point energy of each mode, N{uj)duj is the number of modes with frequencies 
between u and u + du in the presence of boundaries and NQ{uj)duj is the corresponding quantity 
evaluated in empty space. If Ut^{x) are the eigenfunctions of with eigenvalues —uj'^, it is well 
known that the asymptotic distribution of eigenfunctions for large u is given by 

iVH = ^±f^ + -L/ xdS + 0{u-^) (2) 
ZTT 8tt /TT^ JdM 

for Neumann and Dirichlet boundary conditions if the sign is positive or negative, respectively. V 
is the three- volume of M, S is the surface area of dM and x is the trace of the extrinsic curvature. 
Since No(uj) = for scalar fields {E)"^"^^ diverges like uo'^duo. For electromagnetic and 

spinor fields a similar conclusion can be obtained. 
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Baacke and Igarashi investigate the structure of the divergencies of the regularized zero- 
point energy of massive fermions confined in a spherical cavity, assuming the condition of zero 
current across the surface. Considering a massless fermionic field inside and outside a spherical 
shell, Milton showed that there was a cancellation of these divergencies. It is clear that if 
we consider fermions only in the interior of the shell, it is necessary to introduce contact terms 
to deal with the divergencies of the regularized zero-point energy. Thus a question arises: using 
an analytical regularization procedure, is it possible to present a geometric configuration with 
confined fields where the contact terms are not necessary? 

Our purpose is to compute the zero-point energy due to a free massless fermionic field confined 
in the interior of two parallel plates in any number of space-time dimensions. The same idea 
was used by Chodos and Thorn 0, where the Dirac field is confined between two parallel plates 
separated by a distance L in the z direction - a slab. We treat the number of dimensions as a 
continuous parameter and to regularize the zero-point energy we will make use of dimensional 
regularization in the continuous variables (related to the sides of the box with infinite length) and 
then we analytically extend the Hurwitz zeta-function that appears after dimensional regulariza- 
tion. As a consequence, we have been able to prove that the regularized zero-point energy density 
is finite in any number of space-time dimensions. In other words, it is not necessary to introduce 
counterterms in the bare lagrangian, i.e., we have a situation with regularization without renor- 
malization. Some comments are in order. The bag boundary conditions for spin-half fields in 
rectangular cavities, that is, the assumption that the field exists only inside some region, can be 
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implemented only if there is solely one direction of finite extension, because the presence of the 



corners prevents solutions to the massless Dirac equation for higher- dimensional hypercubes ||10|| . 
The situation for massive spinor fields is even worse, in that there is not a solution even in the 
presence of only one direction of finite size |TT]. 



Consider an ultrastatic D-dimensional flat manifold. The zero-point energy of the Dirac field 



is: 

Jd, 



0\H\0) = -c{D) I ^c^„ (3) 



D-2 
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where c{D) is the number of different spin states in a D-dimensional space-time. c{D) = 2 
for D even, and c{D) = for odd D. In the above, 

^l = pi + pl + -+pl-i + pl (4) 

are the eigenfrequencies of the orthonormal modes, basis in the space of solutions of the free Dirac 
equation, i.e., these modes are positive frequency with respect to a timelike Killing vector dt- 

We now turn to the calculation of the zero-point energy of the Dirac field in the D-dimensional 
slab configuration, that is, in the presence of two parallel plates placed in the ci-direction, at 
positions Xd = and Xd = L. The boundary condition that we will impose is that there is no 
particle current through the walls, and consequently we will call this configuration the slab-bag. 
In Lorentz covariant form we have: 

= 0, (5) 
where ri^ = (0, ff), fj being the unit spatial vector normal to the surface and directed to the interior 



of the slab. The boundary conditions are satisfied only if the allowed values for the momentum in 



the d-direction are given by |12 



1 TT 

Pd = Pdin,L) = {n+ -)-, n = 0,1,2,3,... (6) 

The zero-point energy of the Dirac field in the slab-bag configuration, taking into account Eqs. 
(^) and (^), and also noting that the integrations over the momenta in Eq. for the present 
configuration correspond to d — 1 integrations and one summation, is given by: 

(0|^|0) = -c{D) n (^) J dp, J dp, ... J dp,., {pI + ■■■+P>^ . (7) 

j=l Q Q Q n=0 

Defining the total zero-point energy per unit area of the plates, that is, the vacuum energy density: 

Q\H\ 

we arrive at: 



i/ O 

eD= ^^' =ed+i{L), (8) 



eDiL) = --^^,t rd^-^p{pl + ...+plin,L)f\ (9) 

The expression above is clearly divergent, both in the integration and in the summation variables. 
To regularize it we will use dimensional regularization in the continuous variables and then analyt- 
ically extend the Hurwitz zeta-function that will appear after dimensional regularization. Using 
the well-known result of dimensional regularization, i.e.: 

d'^u nif d\ 1 

:^^-n]——ir, (10) 
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it is easy to show that the vacuum energy per unit area is given by: 



(L) = r [-^) Y.{n + (11) 



n=0 



Next we can define f{D) = '^^^^'!^d-i , and note that: 

OO 1 

C(^>?) = Er^lI> 0,-1,-2,-3,... (12) 

n=0 V'' + ^1 

is the Hurwitz zeta-function, which is analytic for Re{z) > 1, in terms of which we can write the 
vacuum energy density of the massless Dirac field in the slab-bag configuration as: 

The two cases, i.e., D even and D odd, are treated separately. We first examine the case of 
even-D. For even-D, the gamma-function presents no poles, and using the principle of analytic 
continuation one can write that: 

a-m,v) = ^^^^, m = 1,2,3,... (14) 
m + 1 

and also the relation between the Bernoulli polynomials and Bernoulli numbers: 5u(|) = —(1 — 
2^^^)Bd, to obtain the vacuum energy density: 

^^(^) = -§Sr(l^) (2^-^-1)^ for D even. (15) 

As a conclusion, Eq. (^) shows us that the regularized zero-point energy per unit area of massless 
fermions in an even- dimensional slab-bag is free of divergencies. In particular, calculating the above 
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expression for D = 4, one finds the value: 



in accordance with 13 . 



In the case where D is odd, one begins with the duphcation formula for the gamma functions: 



n-'^-f^n-s), (17) 



then defines 



2 

in terms of which the vacuum energy density now reads: 



^^(^) = i^r(i-D)c(i-A^). (19) 

Let us now show that the zeta-function times the gamma-function in the equation above is a 
meromorphic function in the whole complex plane. First we write the integral representation of 
the Hurwitz zeta-function: 



oo 



e 



t(l-q) 



Ciz,q) = —-I dtf-'^-, (20) 
1 [z) Jo — 1 

and we separate the range of integration as: [0, 1) + [l,oo), in order to take advantage of the 



expansion: 

- = V :^!ii^r, < |t| < 27r. (21) 
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The integral in the second range, i.e., [l,oo), is a regular function of z, because the integrand 
diverges for small t only, and it will be denoted by 



OO p" 



(22) 



Integrating term by term in the first range, and using Eq. (^2]), we have that: 



Tiz) 



hi{z,q) + J2 



n=0 



:-l)"i^n(g) 1 

n\ z + n — 1 



(23) 



where we have used that -B„(l — g) = (— l)"i?„(g). Gathering everything, we rewrite Eq. (|1^) as: 



9{D) 



ID-I 

where we used again that -Bn(| 



l\ ^ (-l)"(2i-" - 1)5„ 1 



n=0 



n-D 



(24) 



— (1 — 2^ ")-B„. We see from Eq. ( P^ that there is a pole in 
the summation for n = D, with residue: 



ResleoiL)] 



g{D) {-1)^(2''^ -1)Bd 



ID-I 



D\ 



(25) 



But note that for odd-i^, D ^ \, Bo = 0, and the residue of the vacuum energy density in the 
odd-dimensional case vanishes. Although Bi = — |, Eq. ( P^D also states that the residue vanishes 
for D = 1. Eq. ( plf ) is general, and it is valid for the even- dimensional case as well. Although 
there also appears a pole for even-D, it is cancelled by the gamma function in the denominator 
of Eq. ([TBp, and Eq. (|TBp is the result, as can easily be checked. In this way we have proven 
our assertion that the regularized zero-point energy of massless fermions in the slab-bag is free of 
divergencies. 



From Eq. (|1^) one could have taken another route and use that C{s, I) = (2^^ - l)C(s) and the 
reflection formula: 



r[^)as) = 7r^-'/'T[^)ai-s), (26) 



valid for all s, to find that: 



s.W^-S^}^l^r{!l)ao, ,27) 

The appearance of possible divergent terms becomes hidden automatically when one uses the 
reflection formula for zeta functions. The result above is finite for all positive D, and is always 
negative, as shown in Fig. 1. (In the plot of the energy we left aside the spin degeneracy factor 
c{D)). It tends to — oo either for D —>■ and for D oo, taking on the maximum value of 
—4.9 X 10^^ at -D ~ 26.1. For D = 4 result Eq. ([16|) is of course obtained, and this value is 7/2 
times the Casimir energy for a scalar field satisfying Dirichlet boundary conditions on two parallel 
plates. The Casimir energy for the scalar field with Dirichlet boundary conditions as a function of 



the number of dimensions was calculated in Ref . |jTO[ ; it is also always negative and the maximum 
occurs exactly at the same value of D. 

The pressure that the vacuum exerts on the plates is also negative, which means that it acts 
tending to approximate them, and it is given by: 

In conclusion, in this letter we discuss open questions concerning the divergent pieces of the 
regularized zero-point energy of a fermionic field confined in a finite volume. We show that the 
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regularized zero-point energy of massless fermions in the interior of two parallel plates, the slab- 
bag, is finite for any number of space-time dimensions. For such, we made use of dimensional 
regularization in the continuous momenta, related with the d — 1 transversal directions of infinite 
length of the bag, and then we analytically extended the Hurwitz zeta-function thereby obtained. 
For even-dimensional space-times the zero-point energy is finite, without need of renormalization; 
in the odd-D case it is shown that the residue of the polar part vanishes identically and, therefore, 
the usual subtraction of the polar part of a divergent quantity is also not necessary. Also using 



an analytic continuation procedure, Dolan and Nash |jT^ compute the Casimir energy of massless 
spin zero fields on spheres S*^ and find that, for the case of conformal coupling, the Casimir 
energy vanishes for even-dimensional spheres, presenting no divergencies, and for odd- dimensional 
spheres the residue of the pole which appears is shown to vanish giving a finite value for the Casimir 
energy. Self-interacting fermions at finite temperature and density was studied by Afianos et al. 
Tsj . These authors investigate, at the (large A^) approximation, the behavior of the effective 



m 



potential of the Gross-Neveu model. Still using a mix between zeta and dimensional regularization 
procedures, they found that the regularized effective potential is finite for odd-dimensional space- 
times. In the case of the Yukawa model, the same kind of behavior of the regularized effective 
potential was found [|^. 

One might be tempted to suppose that these results are related with the work of Svaiter and 



Svaiter [|1^. These authors proved that the zeta function method is equivalent to the cut-off 
method with the subtraction of the polar terms for a neutral massless scalar field. The basic 
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idea employed by Casimir |T|], who introduced the cut-off method, is that although the zero-point 
energy is divergent, changes in the configuration lead to a finite shift in the total energy. In the 
exponential cut-off method, in order to evaluate this shift, the total energy is regularized before the 
subtraction from the energy of a reference configuration. This total energy is obtained by adding 
the regularized zero-point energy of the field inside and outside the cavity. This approach seems 
very natural if we are dealing with a system in which there is field inside and outside the cavity. 
However, in the bag configuration, where the field is supposed to exist only inside some region, this 



is not a natural approach and we follow the conclusion of the authors of Ref. IjT^ that the Casimir 
energy evaluated using the analytic extension of the zeta function throws away these divergent 
terms. A natural extension of this work is to calculate the regularized zero-point energy for the 
fermion field for different confining geometries, such as the cylinder and the sphere, still using an 
analytical regularization procedure. Now, whether or not the analytical regularization procedure 
throws away the polar part of the regularized vacuum energy deserves further investigation. 
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